/****************************************
Table A.25
*****************************************/

* Load Census Data Cleaned and Collapsed
use "$dir/Data/Final/final_collapsed_clean.dta", clear
	
gen i = birthstatefip
gen t = (birthdecade/10) + 1 // someone who is 10 when ban is put in place is treated (time 0)
gen Ei = floor(year_of_ban/10) 	// decade when unit is first treated
gen K = t-Ei 			// "relative time" num decades since treated (could be missing if never-treated)

gen lastcohort = (Ei >= 193 - 1 ) & Ei != . /* ban after 1920 */



forvalues l = 1/6 {
	gen A`l'event = K==-`l'
}

forvalues l = 0/6 {
	gen P`l'event = K==`l'
}

gen insample6treat = abs(K) <= 6 & K != .



* Interaction treatment vars
forvalues pre = 6(-1)2 {
	cap n gen cmH_Pre`pre' = cmH * A`pre'event
}

forvalues post = 0/6 {
	cap n gen cmH_Post`post' = cmH * P`post'event
}

forvalues x = 0/6 {
	cap label var cmH_Pre`x' "Pre `x'"
	cap label var cmH_Post`x' "Post `x'"
}



drop if year_of_ban > 1940
drop if year >= 1930 /* Drop time-periods when last-treated units have received treatment */

keep if insample6treat == 1

**************************************************
** Transforming and labelling outcome variables **
**************************************************

* Transformations
gen ln_citypop = ln(citypop)
gen ln_occscore = ln(occscore)
gen ln_inchat = log(inchat_levels)
drop inchat_levels
gen occ_disp = 1 - occ_disp_hhi 
gen occ_disp_10digit = 1 - occ_disp_hhi_10digit 
gen occ_disp_occ_class = 1 - occ_disp_hhi_occ_class
gen geog_disp_state = 1 - geog_disp_hhi_state
gen geog_disp_county = 1 - geog_disp_hhi_county
gen geog_disp_enumdist = 1 - geog_disp_hhi_enumdist
gen share_ch_uncmn = 1 - share_ch_cmn 


* Scaling genetic outcomes
foreach var of varlist(gq_inst-genetic_d){
	qui: replace `var' = `var'*10000
}

* Labels
label var birthdecade "Birth decade"
label var cmH ">10% cousin marriage rate (-1858)"
label var birthstatefip "Birth state"
label var resstatefip "Residence state"
label var year "Census year"
label var count_dor "Num census individuals in cell"
label var occscore "Occscore"
label var inchat "LIDO Income (Original)"
label var ln_inchat "LIDO Income"
label var citypop "City size"
label var urban "Urban"
label var farm "Residence: Farm"
label var mig_life "Interstate Migration"
label var mig_life_noban "Interstate Migration (State without ban)"
label var mig_life_ban "Interstate Migration (State with ban)"
label var ln_citypop "Urbanization (log residence pop'n)"
label var divorced "Divorced"
label var mult_genh "Multigenerational HH"
label var nchild "No. of resident children"
label var nchlt5 "No. of resident children (under 5)"
label var agemarr "Age of Marriage"
label var durmarr "Duration of Current Marriage"
label var ncouples "No. of couples in unit"
label var ncouples_wo "No. of couples in unit (w/o own)"
label var nfams "No. of unrelated families in unit"
label var famsize "Family Size"
label var nsibs "No. of siblings in unit"
label var nm_single "Never Married/Single"
label var occ_disp_hhi "Occupational Dispersion HHI"
label var occ_disp_hhi_10digit "Occupational Dispersion HHI (10-digit codes)"
label var occ_disp_hhi_occ_class "Occupational Dispersion HHI (Farmer, White Collar, Blue Collar, Self-Employed)"
label var prank "Income percentile-rank"
label var geog_disp_hhi_state "Geographical Dispersion (State)"
label var geog_disp_hhi_county "Geographical Dispersion (County)"
label var geog_disp_hhi_enumdist "Geographical Dispersion (Enumeration District)"
label var occ_disp "Occupational Dispersion (HHI)"
label var occ_disp_10digit "Occupational Dispersion (HHI) (10-digit)"
label var occ_disp_occ_class "Occupational Dispersion (HHI) - White/Blue Collar, Farmer, Self-Employed"
label var geog_disp_state "Geographical Dispersion (HHI) - State"
label var geog_disp_county "Geographical Dispersion (HHI) - County"
label var geog_disp_enumdist "Geographical Dispersion (HHI) - Enumeration District"
label var gq_inst "Institutionalized, per 10,000"
label var gq_inst_med "Institutionalized (Medical), per 10,000"
label var gq_poorhouse "Institutionalized (Poorhouse), per 10,000"
label var hmem_disable "Household member (Disabled), per 10,000"
label var blind_deaf "Blind/Deaf, per 10,000"
label var genetic_d "Institutionalized (Medical) + Disability, per 10,000"
label var share_ch_cmn "Share childr. with common names"
label var share_ch_uncmn "Share childr. with uncommon names"

foreach v of varlist(ln_occscore_nm-genetic_d_nm){
	label var `v' "Number of individuals with non-missing outcome in cell"
}

compress

/**************************************
Table A.25
**************************************/

eststo clear
eventstudyinteract ln_citypop ///
	cmH ///
    cmH_P* ///
	[aweight = ln_citypop_nm] , absorb(i.i##i.t i.year  i.i##i.cmH) vce(cluster i) ///
	cohort(Ei) control_cohort(lastcohort)
	
 	matrix b_iw = e(b_iw) // Save the results for table output
 	matrix V_iw = e(V_iw)
	capture matrix V_iw = diag(e(V_iw))
 	erepost b = b_iw V = V_iw // ssc install erepost
 	esttab
	eststo A1
    sum ln_citypop [fw = ln_citypop_nm]
	local obs1 = r(N)
    su ln_citypop [aw = ln_citypop_nm] if e(sample) == 1
    local meanA = string(round(r(mean), .01))
 
 
eventstudyinteract ln_occscore ///
	cmH ///
    cmH_P* ///
	[aweight = ln_occscore_nm] , absorb(i.i##i.t i.year  i.i##i.cmH) vce(cluster i) ///
	cohort(Ei) control_cohort(lastcohort)
	
 	matrix b_iw = e(b_iw) // Save the results for table output
 	matrix V_iw = e(V_iw)
	capture matrix V_iw = diag(e(V_iw))
 	erepost b = b_iw V = V_iw // ssc install erepost
 	esttab
	eststo A2
	sum ln_occscore [fw = ln_occscore_nm]
	local obs2 = r(N)
    su ln_occscore [aw = ln_occscore_nm] if e(sample) == 1
    local meanB = string(round(r(mean), .01))
 

 estout A1 A2 using "$dir/Output/TablesFigures/results_lasttreated.tex", style(tex) replace ///
    keep(cmH_*) ///
	cells(b(star fmt(%9.4f)) se(par)) ///
	nolabel collabels(none) mlabels(none) starlevels(* 0.10 ** 0.05 *** 0.01) ///
	varlabels(cmH_Pre6 "-6" cmH_Pre5 "-5" ///
	cmH_Pre4 "-4" cmH_Pre3 "-3" ///
	cmH_Pre2 "-2" cmH_Post0 "0" ///
	cmH_Post1 "1" cmH_Post2 "2" ///
	cmH_Post3 "3" cmH_Post4 "4" ///
	cmH_Post5 "5" cmH_Post6 "6")

local tex " \\ \hline"
local tex "`tex' Observations & `obs1' & `obs2' \\"
local tex "`tex' Mean Dep. Var & `meanA' & `meanB'   \\"
local tex "`tex' State X Birth Decade FE & Yes & Yes \\"
local tex "`tex' State X High CM F.E. & Yes & Yes  \\"
local tex "`tex' Census year F.E. & Yes & Yes  \\"
local tex "`tex' High CM X Birth Decade Trends & Yes & Yes  \\"
local tex "`tex' \multicolumn{3}{p{8cm}}{\tiny \textit{Notes:} Standard errors clustered at state-level."
local tex "`tex' *** p<0.01, ** p<0.05, * p<0.1.} \\ \end{tabular} }"
	
esttab A1 A2 using "$dir/Output/TablesFigures/results_lasttreateda.tex", style(tex) replace booktabs ///
	d(*) nolabel collabels(none) noobs postfoot("`tex'") nonum ///
	mtitles("(1)" "(2)") ///
	mgroups("Urbanization (log residence pop'n')" "Occupational Income (Log)", pattern(1 1) ///
	prefix(\multicolumn{@span}{c}{) suffix(}) span erepeat(\cmidrule(lr){@span}))
	

